Microscopic origin of low frequency flux noise in Josephson circuits 
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We analyze the data and discuss their imphcations for the microscopic origin of the low frequency 
flux noise in superconducting circuits. We argue that this noise is produced by spins at superconduc- 
tor insulator boundary whose dynamics is due to RKKY interaction. We show that this mechanism 
explains size independence of the noise, different frequency dependences of the spectra reported in 
large and small SQUIDs and gives the correct intensity for realistic parameters. 
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Studies of the flux noise in superconducting structures 
have a long history that began in 80s with the demonstra- 
tion that it is the flux and not the critical current noise 
that limits the sensitivity of dc SQUIDs (Superconduct- 
ing QUantum Interference Devices) [l|, Q . The noise phe- 
nomenological characterization performed at this time re- 
vealed many puzzling features that defled a simple model 
of its microscopic origin, so this problem was put aside 
and largely forgotten. Recently, the interest to it was re- 
newed when it was realized that the flux noise limits the 
coherence of qubits based on superconducting devices Q 
and that the dephasing of 'flux' qubits is due by low fre- 
quency flux noise with intensity comparable to the one 
measured in dc-SQUIDs 0, H] ■ It is likely that the same 
flux noise will limit the quantum coherence in 'phase' 
qubits 0. 

Recently two models for the excess 1// flux noise were 
proposed. The first onef?! proposes that the flux noise is 
due to the electrons that hop between traps in which their 
spins have fixed, random orientations. The second model 
[8| attributes the noise to the electrons that experience 
spin fiips induced by the interaction with tunneling Two 
Level Systems (TLSs) and phonons. Both models rely on 
some assumptions that are difficult to justify: for exam- 
ple, in order to match the intensity of the noise spectral 
density reported in the experiments, the number of ther- 
mally activated TLSs present in the oxide layer has to be 
much larger than in a typical glass at the same temper- 
ature. For instance, in a typical loop of radius R — Ifim 
and volume lO^nm'^ the observed noise value implies ac- 
tivation of 10^ — 10^ spin fluctuators with magnetic mo- 
ment fiB while a typical glass of the volume has about 
10 thermally excited TLSs at T = O.IK. This Letter 
has two goals: (i) to present a critical analysis of the flux 
noise phenomenology and its implication for the possible 
models of its microscopic origin, (ii) to propose a novel 
mechanism in which the low frequency noise is due to 
spin diffusion on the superconductor surface generated by 
the exchange mediated by the conduction electrons. We 
demonstrate that this spin dynamics together with the 
spatial dependence of the surface current density on the 
thin superconducting SQUID loop leads to low frequency 



1//" flux noise spectral density with a G {0, 1} and that 
the intensity of the noise does not depend on the area of 
the SQUID, as long the ratio R/W remains constant {W 
denotes the width of the SQUID line); however, the de- 
tails of diffusion in large SQUIDs, with W ^ lOO/im and 
small SQUIDs having W ~ Ijim might be different. In 
particular, the frequency dependence of the noise spec- 
trum might vary depending on the size of the SQUIDs 
and the measured frequency range. 

All experiments agree on the magnitude of the 
noise at frequency / ~ IHz and its area indepen- 
dence. Specifically, Wellstood observed noise spec- 
tra 5y^(lHz) « (4- 10)Ai*oHz"^/^ at temperature be- 
low O.IK in Nb and Pb dc SQUIDs with sizes in the 
range R,W ^ (30 — 300)^m on Si/SiO substrate with 
Nb/NbO/PnIn Josephson junctions. Cromar et al Q 
reported the value ^^^(IHz) = 2.3^$oHz-^/2 for Nb 
SQUIDs device with very high quality Nb/AlO/Nb junc- 
tions at 4K. Finally, Bialczak et al [l^ measured 
5y^(lHz) = 2ii^oHz-^l^ at 20mK in a Al SQUID loop 
of size W ~ l/.tm on sapphire substrate with Al/AlO/Al 
Josephson junctions. Although the details of temper- 
ature dependence were studied only in 0, all experi- 
ments agree that the noise does not decrease at very 
low temperatures. Similarly, all data show homoge- 
neous noise spectra where single strong ffuctuators can- 
not be resolved. The frequency dependence of the noise is 
more controversial. Namely, Wellstood's flux noise power 
spectra in the frequency range (1 — 10'^)Hz displayed 
1//" dependence with exponent a = 0.66 at low tem- 
peratures 0.022K < T < IK and with exponent a = 1 at 
IK < r < 4.2K. The data Q show the dependence f-° '^ 
in the interval (400- 10^) Hz at 4K, below 400 Hz the 
frequency dependence decreases to approximately /~°'^ 
and completely ceases in the (0.1 — 40) Hz interval. Bial- 
czak et al reported 1//" spectrum with a = 0.95 in 
the frequency range (10~^ — l)Hz at 20mK. The tem- 
perature dependence of the noise 0] shows two different 
temperature regimes: at < 0.5K the noise is T indepen- 
dent, while at IK < T < 4.2K it displays dependence 
with the crossover regime (0.5K < T < IK) that is non- 
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monotonous in some samples. Two distinct regimes sug- 
gests two different microscopic mechanisms for the noise 
at low and high temperature; in the following we shall 
focus on the low temperature regime. As it will be clear 
below, a very important piece of information is the high 
frequency cutoff of the 1//" dependence. Unfortunately, 
no direct measurements are available but the observed 
dephasing of the flux qubits indicates that this cutoff is 
at least 10MHz[il|. 

We now discuss the implications of the data for the 
noise origin. The noise persistence at low temperatures 
indicates that it is due to a subsystem characterized by 
very low energy scales, smaller than the minimal temper- 
ature available experimentally {'^ 20mK) . This rules out 
the thermally excited TLSs [7|] or vortices and points to- 
wards weakly interacting nuclear or electron spins. The 
spin mechanisms agrees also with the observation of ho- 
mogeneous low frequency noise power spectra (which is 
incompatible with the vortex origin). Nuclear spins can 
be excluded for three reasons: first, the flux produced 
by each spin scales as l/L (where L ^ R,W is the lin- 
ear size of the device) thereby leading to Slf^{Lu) oc L^^'^ 
while data are roughly size independent; second, all fre- 
quency scales associated with nuclear spins are very low 
(/ < IkHz) in contrast with the results of the dephas- 
ing analysis which shows that 1// persists up to lOMHz 
ll|; third, one expects that nuclear spin noise would be 
substrate dependent 0]. Paramagnetic electron spins lo- 
cated on the superconductor or insulator interfaces seem 
to be more promising candidates since their contribution 
to the flux noise is roughly size independent. Properties 
of these spins were extensively studied for Si/Si02 inter- 
faces. ESR experiments have shown that: (i) the surface 
density of spins varies between 1^20 ~ 10"^° — lO^^cm"^ 
[1^ : (ii) the g factor of these spins is isotropic and it 
has value g = 2.00136 ± 0.00003 [i3|. As we show below, 
such surface density is barely sufficient to explain the 
level of flux noise at IHz if one assumes that all these 
spins remain active at low temperatures but is difficult to 
reconcile them with schemes in which only a small per- 
centage of the spins remain active at low T (3] ■ The value 
g = 2 shows that the spin orbit coupling is very weak in- 
dicating that the interaction between paramagnetic spins 
and TLSs is very small in contrast to assumptions of Rcfs. 

The dynamics of the electron spins in the insulator 
substrate is due to the interaction with other electron 
spins or with surrounding nuclei. For a dilute spin 
system such as Si/Si02 interfaces all energy scales as- 
sociated with these interactions are too small to ac- 
count for the wide frequency range observed experimen- 
tally: both dipole-dipole interaction between electron 
spins with density U2d ~ lO^^cm"^ and their interaction 
with nuclear moments of Si that have natural concen- 
tration of 5% correspond to / w lOkHz. Estimating the 
total flux noise produced by these spins, i.e. J S{Lu)dLL} we 



get - i^lo^^B) {R/W)iy2D^l with (AtoMs) ~ IQ-^^cm^ 
which is of the right order of magnitude but somewhat 
smaller than the observed noise. We conclude that these 
spins in the insulator arc unlikely to provide the domi- 
nant source of noise. 

The energy scales are much larger for the electron 
spins in the proximity of the superconductor which al- 
lows RKKY interaction between them which is a much 
stronger coupling than dipole-dipole interaction in the in- 
sulator. This interaction is due to the virtual scattering 
of conduction electrons off a magnetic impurity described 
by the Kondo Hamiltonian: Hk = J^S ■ a where S is the 
spin operator for the impurity, a is the spin operator of 
a conduction electron and S is the exchange constant. 
Integrating out the conduction electrons one gets: 



RKKY 



(1) 



where V{r)~VQ{r)e "^^^^r"^ cos^p and changes 
quickly on the length scale of the Fermi wavelength \p 
[l4|. The interaction 'constant' Vo(r') is a weak function 
of the distance, it is controlled by the electron density of 
states V and Kondo temperature: Vo(?') = (27r)~^z/J'^(7'), 
with J{r) = 2[v\t? [vp / {tTk))]'^ so that the average in- 
teraction at r <C ^ 



{V\r)f" = 



(2) 



Because of this interaction, the magnetization M{t,r) 
of spins averaged over the volume that contains ^ 1 
spins obeys the diffusion equation: 



d 
dt 



M{t, r)=0 



(3) 



with diffusion coefficient V which depends on the typ- 
ical distance between the spins on the surface, i.e 
= \Jv2B ~ (10 — 10^) nm and the average interaction 

(y^(r))^^^ CJ- Typical electron density of states for 
Al, Pb and Nb are respectively: vx\ = 35/eFnm'^, 
z/pb = 44/eT^nm'^ and j^Nb — 160/el/nm'^. Assuming 
Kondo temperatures w 0.01 — IK, we estimate: 
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This model neglects a few important physical effects. 
First, it neglects the spin orbit scattering and assumes 
that the diffusion process involves only electron spins lo- 
cated on the SI interface. As a result, the total mag- 
netization M of the spins in contact with the supercon- 
ductor is conserved. Second, the estimate for V Eq.Q 
assumes that the spins are in direct contact with the 
metal. However, paramagnetic spins responsible for the 
flux noise are likely to be located in the surface oxide of 
the thickness c? = (2 — 3)nm with some of them further 
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away from the superconducting wire. For impurity lo- 
cated at depth y from the superconductor the strength 



of RKKY interaction decreases as V{r, y) 



V{r), 



where is the atomic distance. A more reahstic model 
should include 'fast' spins at the surface with the diffu- 
sion constant given in Eq. (g]) and slower spins coupled 
to the 'fast' subsystem by a weakened RKKY interaction. 
Third, the diffusion approximation for the spin dynam- 
ics neglects the effect of the rare pairs of spins located 
at distances much smaller than the average distance be- 
tween the spins. Such spins are strongly coupled with 
each other, the difference in the energy of their triplet 
and singlet state is much larger than their coupling to 
their neighbors, so they change their state rarely. This 
mechanism generates an additional noise at low frequen- 
cies. 

In order to find the effective flux $e// produced by 
the spin magnetization we determine the spin energy E 
in the field of the test current / in the loop. We find that: 



dE 
IT 



(5) 



Here is the Bohr magneton, S{r) is the spin density 
operator and B(r) denotes the probing magnetic field. 
Conservation of the total magnetization by spin diffusion 
means that it would not produce any noise if the probing 
magnetic field were uniform. In fact, it is not: the SQUID 
loop is typically a strip conductor of width W greater 
than its thickness with length L ^ W. For the film 
thickness less or comparable with the penetration depth 
A, the dependence of the current density on x near the 
center of the strip is J,(x) = 2I/{TrW)[l - {2x/Wf]-^/^ 
for — W/2 + X < X < W/2 + A, while the current de nsity 
falls away exponentially to zero at the edges ±W/2 [15|. 
This current density results in a probing magnetic field 
B{x) ~ ^Js{x). The spin diffusion together with the di- 
vergency of the surface current density close to the edges 
of the loop generates 1// flux spectral density: 



l-w/2 

($^$o) = ig^J'Bf L I dxdx' 

J -w/2 



,Sr{x)B{x)So{x')B{x') 



p 



(6) 

where x is the coordinate across the wire strip. To com- 
pute the integral in Eq. ([6]) we expand the spin density 
operator St (x) as a series of orthonormal eigenfunctions 
of the diffusion equation ([3]) with boundary conditions 
ensuring zero magnetization current at the wire bound- 
ary, i.e. ^Mtix = ±W/2) = 0: 



dt 



St{x) 



— ^^Wcos 



W 



(7) 



q—7Tn/W 

where n is positive integer. By substituting the expansion 



given in Eq. ^ into Eq. ([6]) we find: 

($.$0) = ig^iBrLY,B^^{S,iT)S,m 



(8) 



where 
B„ 



w/2 



B{x) 
dx — - — cos 



w/2 ^ 

/^O ^ (qW\ fqW 



w 
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where j7o(a;) is the Bcssel function. The Fourier trans- 
form of the spin density correlator ([8|) is found from the 
solution of the diffusion equation ([3|): 



2 ^2 + (Vq^)^ 



(10) 



where as ~ pd is the surface density of paramagnetic 
spins ^. We can define two frequency regimes: small 
frequencies with / ^ fw and large frequencies with 
/ ^ fw- fw is the characteristic equilibrium frequency 
for spins that diffuse across SQUID of width W: 



w 



V 



10-2 - 10-iHz if VF - 100/im; 
10^ - lO^Hz if VF - 1/im. 



At small frequencies, the flux noise spectrum given in Eq. 
^ is white, with noise amplitude given by: 



(11) 



reduces to B„ 



\ 2tt J 'W fw 
: large frequen( 
and the flux noise spectrum 



where Jo{Tr) = -0.3042. At large frequencies, Eq. ^ 
^/io 1 



ttW ^ 



given in Eq.jS]) becomes: 
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(MoMb) 



L 
R 1 



dq Vq^ 
~Tuj^ + {Vq^y 



(12) 



where we have written explicitly the length of the SQUID 
loop L = 2ttR. At intermediate frequencies, we expect 
a crossover between 1// and white noise behavior. It is 
quite straightforward to estimate the intensity of the 1// 
noise. Assuming that (Ts « lO^^m^^ (similar density was 
reported for Si/SiOj interfaces) and that R/W ^ 10 we 
flnd flux noise spectral density S'$(lIIz) w 3(/i<f>o)^IIz~"^j 
in agreement with the observed "universal" value. Be- 
cause the noise is due to the spins on the surface, its 
level has the same R/W dependence as in Ref [lo| . 

Thus, the spin diffusion model explains the excess flux 
noise measured in large SQUIDs QjSIj which spectra 
correspond to the intermediate/high frequencies regime. 
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but not the 1// noise observed in much smaller de- 
vices since the latter was measured in the range 
corresponding to the low frequency regime where purely 
spin diffusion model predicts a constant spectral density. 
However, two physical effects missing in the model that 
were mentioned above are very likely to produce a sig- 
nificant low frequency noise in the smaller SQUIDs: the 
presence of weakly coupled spins further away from su- 
perconductor and the presence of strongly coupled spin 
pairs. Indeed, assuming flat distribution of the spin 
depth inside the insulating layer one gets an exponential 
distribution of the coupling to the spins on SI interface 
P{J) oc 1/J that directly translates into the 1// spec- 
trum of the noise generated by these spins. The intensity 
of this noise is determined by the areal density, V2JJ of the 
spins in the layer of approximately atomic depth ~ 2ao. 
Generally, one expects ~ ^20] a significant difference 
in the values of j/j^j, V2D might lead to a complicated fre- 
quency dependence of the noise observed in a wide fre- 
quency range: e.g. a 1// behavior at high frequencies 
(due to 'fast' spin diffusion) followed by a partial satura- 
tion at intermediate frequencies which is followed again 
by 1// regime at very low frequencies due to deep and 
weakly coupled spins. Our preliminary analysis shows 
that close pairs of spins strongly coupled to each other 
by RKKY interaction also lead to 1// contribution to 
the low frequency noise. 

Finally, we discuss experimental tests of the proposed 
model. The crucial ingredient of our analysis is the rough 
temperature independence of the noise below 200 mK [2] , 
it would be important to verify it for small devices. The 
spin origin of the noise can be tested by applying a signif- 
icant external magnetic field. If this field is larger than 
the local field, i?ioc, produced by the spin neighbors, the 
spin rotates around the axis determined by the external 
field. If it is orthogonal to the probing field, fast rotation 
of the spin implies that the effective spin noise is shifted 
to high frequencies. If these fields are parallel, the effect 
is much less. The effective probing field acting on the 
spins on the insulator boundary inside the SQUID loop is 
mostly perpendicular to the surface of the sample, while 
the probing field of the spins on SI boundary is parallel 
to it. Thus, applying magnetic field in different direc- 
tions one can verify the spin mechanisms and determine 
the spin location. The local field that should be exceeded 
in these experiments is of the order of _Bioc ^ lOOG for 
the spins lOnm apart on SI surface and of the order of 
B\oc ^ O.IG for the spins in the insulator. Random po- 
sition of spins in these models implies that there will be 
always strongly coupled spin pairs capable of producing 
the low frequency noise but the number of such pairs 
should go down rapidly with field. The validity of the 
spin models discussed in this paper can be also tested di- 
rectly by fabrication of the samples with decreased den- 
sity of spin defects on the surface of the insulator and 
by protecting the surface of superconductor of a layer of 



another metal, e.g. Re. In this paper we have not discuss 
complicated mechanisms involving the combined effects 
of electron and nuclear spins such as electron spin rota- 
tion induced by nuclear spin nearby. We believe that it is 
unlikely that these mechanisms can produce sufficiently 
high upper frequency cutoff and sufficient noise level for 
a natural Si with a low concentration of nuclear spins 
but this should be also verified experimentally by mea- 
suring noise on isotopically pure Si substrates. Finally, 
in our model the 1// dependence of the noise is due to 
'fast' diffusion and the divergent dependence of the probe 
magnetic field at the edge of the SQUID loop. Thicker 
SQUID loops have different spatial current distribution, 
this should affect the frequency dependence of the noise. 
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